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Quantifier 



Data members 

(required) 

Array dimensions /, /, K 
Scalar type 

Density array (address) 
IVIethod-choice flags 

(as needed) 
Bound_density array 
Free_density array 
Binding_rate array 
Creation^rate array 
/x/x/T segmentation array (address) 
Segment-surface list (address) 
Transport model (address) 



635 



Function members 

(scalar-valued) 

VolumeJn_range( [Co], [Ci], [segment] ) 
lntegrate( [segment | triangulated_surface] ) 
Efflux{ triangulated_surface ) 

(array-valued) 

Difference( subtrahend, [segment], [divisor], [target] ) 
Multijmage( arrayjist, algorithm, [target] ) 
lmage_process( operator, [segment], [target] ) 



+ member-setting functions 

J 
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TemporalQuantifier 



711 ] Data members 
(required) 

^Array dimensions /, J, K 
^Number iV of times 
Times /i, tz, , tpi 
Scalar type 
Density array list 
Method-choice flags 

(as needed) 

iound density array list 
-Free density array list 
binding rate array list 
-Creation rate array list 
-/x/xA^ segmentation array (static or list) 
-Segment-surface (static or list) 
^Transport model (address) 
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Function members 

(scalar-valued) 
^VolumeJn_range( [Co], [Ci], [segment], [times]) 
'lntegrate( [segment | triangulated_surface], [times] ) 
'Efflux( triangulated_surfacejist, [times] ) 

(array-valued) 

Rate_of_change( [segment], [divisor], [target], [times] ) 
-Multijmage( arrayjist, algorithm, [target], [times] ) 
Jmage_process( operator, [segment], [target], [times] ) 



73/ 



+ member-setting functions 
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Initialize pointer on 
first triangle. 
Volume6 = 0.0 



810 



Volume += prism volume 



last triangle? 



820/-\_y 
821v^ 



No 



next 
triangle 



Yes 



850 



Return voiume6/ 
6.0 



910 

911 
912 

920 
921 
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Initialize pointer on 
first gridpoint. 
Count = 0 . 



Ci< c(point)< C2? 



Yes 



In segment S? 



Yes 
T 



No 



Count += 1 



No 



Last grid point? 



No 



next 
point 



Yes 



Return count xH 
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1040 
104 



1050 



Initialize pointer on first 
volume element. 

V = 0.0 



Find 8-corner type T with respect to 
above/below threshold C 



Interpolate edge lengths relevant to T 



V += polynomial relevant to T 




Last element? 



No 



next 
element 



Yes 




FV 



1010 



1020 



1030 



1140 
1141 



Initialize 'perhaps In S' 
volume element iterator. 

V = 0.0 



FV 



All corners outside segment S? 



All corners Inside S? 




Yes 



Yes 
I 



V += fractional content 



Find 8-corner type T with respect to 
condition *above/below threshold C 



Interpolate edge lengths relevant to T 



V += polynomial relevant to T 



.1110 



1120 



.1130 



Last element? 



No 

T 



next 
element 



Yes 



Return VxH 



1150 



Input array parameters, 

range limits Ci, C2, 
optionally segment S. 



Set range 
condition 

<C2 



FV or FV 



1 

Set r 
cond 
< 


ange 
ition 

C2 






FV or FV 




Return difference 
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1250 



Estimate dc/dt over grid, 
in array with address D 



Select a boundary B where 
c is negligible 



Fix array for free-density values, 
with address A. 



Pass D, A and B to transport 
model inversion routine. 



Return A 
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Create array e, 
initialize values to 0.0 



Triangle index 1 = 0 



1410 



Compute normal vector n of triangle i 



Compute mid-point w of triangle i 



Pass w to transport model 



Receive vector F from transport model 



e[N] += n. F, e[i] = n • F/(triangle area) 
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Last triangle? 



No 



index i += 1 



Yes 



Return A 



1420 



1430 



1440 



1441 



1450 
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